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Abstract 

The interplay between Rashba, Dresselhaus and Zeeman interactions in a quantum well sub- 
mitted to an external magnetic field is studied by means of an accurate analytical solution of the 
Hamiltonian, including electron-electron interactions in a sum rule approach. This solution allows 
to discuss the influence of the spin-orbit coupling on some relevant quantities that have been mea- 
sured in inelastic light scattering and electron-spin resonance experiments on quantum wells. In 
particular, we have evaluated the spin-orbit contribution to the spin splitting of the Landau levels 
and to the splitting of charge- and spin-density excitations. We also discuss how the spin-orbit 
effects change if the applied magnetic field is tilted with respect to the direction perpendicular to 
the quantum well. 

PACS numbers: 73.21.Fg, 73.22.Dj, 73.22.Lp 
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I. INTRODUCTION 



The study of spin-orbit (SO) effects in semiconductor nanostructures has been the object 
of many experimental and theoretical investigations in the last few years, see e.g. Refs. 1-16 
and Refs. therein. In spite of this, the extraction from measurements of the effective spin- 
orbit coupling constant of both Dresselhaus 17 and Bychkov-Rashba 18 ' 19 SO interactions is not 
a simple matter, since the SO corrections to the electron energy spectrum in a magnetic field 
(B) are vanishingly small because they correspond to second order effects in perturbation 
theory. Thus, few physical observables are sensitive enough to the SO interactions and allow 
for a quantitative estimate of their coupling constants. One such observable is the splitting of 
the cyclotron resonance (CR), which has been determined in transmission experiments with 
far-infrared radiation, 20 and is due to the coupling between charge- density and spin-density 
excitations. 21 A less clear example is the change in the Larmor frequency -spin splitting. 22 
The spin splitting has been observed in electron-spin resonance 23 ' 24 and in inelastic light 
scattering experiments. 25 ' 26 

In this work we extend our previous results 21,22 by obtaining an approximate, yet very 
accurate, analytical solution of the quantum well SO Hamiltonian that contains both Dres- 
selhaus and Bychkov-Rashba interactions. In the limit of high magnetic field, this solution 
coincides with the results of second order perturbation theory, and allows to study the SO 
corrections to the Landau levels in a simple way, and to study the transitions induced by an 
external electromagnetic field acting upon the system. 

This work is organized as follows. In Sec. II we present the general formalism for the 
single-particle (sp) Hamiltonian. These results are used in Sec. Ill to study the transitions 
caused by an external electromagnetic field. The role of the electron-electron (e-e) interaction 
is discussed In Sec. IV within a sum rule approach. In Sec. V we discuss the splitting of 
the Landau levels and the appearance of charge- and spin-density modes making, whenever 
possible, qualitative comparisons with the experimental results. 20 ' 24 ' 27 ~ 29 A brief summary is 
presented in Sec. VI, and the generalization of some of the expressions derived in Sec. II to 
the case of tilted magnetic fields is presented in the Appendix. 
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II. SINGLE-PARTICLE STATES 



In the effective mass, dielectric constant approximation, the quantum well Hamiltonian 
H can be written as H = H + ^ X^ J=1 | r . j^q , where if is the one-body Hamiltonian 
consisting of the kinetic, Zeeman, Rashba and Dresselhaus terms 

N N 



i=i i=i 



p+ p - + p-p+ x 



+ -g*fj, B Ba z 



Am T 

(i) 



m = m*m e is the effective electron mass in units of the bare electron mass m e , P ± = P x ±iP y , 
a± — a x ± ia y , where the cr's are the Pauli matrices, and P = —iTiV + represents the 
canonical momentum in terms of the vector potential A, which in the following we write in 
the Landau gauge, A = B(0, x, 0), with B = V x A = Bz. The second term in Eq. (1) 
is the Zeeman energy, where hb = he/(2m e c) is the Bohr magneton, and g* is the effective 
gyromagnetic factor. The third and fourth terms are the usual Rashba and Dresselhaus 
interactions, respectively. Note that for bulk GaAs, taken here as an example, g* = —0.44, 
m* = 0.067, and the dielectric constant is e = 12.4. To simplify the expressions, in the 
following we shall use effective atomic units ft = e 2 /e — m — 1. 
Introducing the operators 

a± = -^=P ± (2) 



with [a~,a + ] = 1 and u c = eB/c being the cyclotron frequency, the sp Hamiltonian h can 
be rewritten as 

h /uj c = \{a + a~ + a~a + ) - -— o z - ]-i~X R (a + a_ - a~a + ) + - X D (a + a + + a~<7_) , (3) 
Z Z 00c Z Z 



where ojl = \g*[J>BB\ is the Larmor frequency and X RD = Xr^^J-^-. For the spinor \<p) = 

( 02 ) ^ We sna ^ use ^ or ^ e com P onen t an( i '2' for the bottom component of any 
spinor), the Schrodinger equation ho\<f>) — e\<f>) adopts the form 



= . (4) 



|(a + a + a a + ) — uj l /(2uj c ) — e iX^a + A^a + 

—iX R a + + Xr>a~ |(a + a~ + a~a + ) + ujl/{2uj c ) — e 

We expand 0i and 2 into oscillator states \n) as 4>i — J2^=o a n\n), 02 = zZ^=o^n\ n ) , 
on which a + and a~ act in the usual way, i.e., |(a + a~ + a~a + )\n) — (n + \)\n), a + \n) 
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), 02 = 



\Jn + l|n+l), a \n) = y/n\n — l), and a |0) = 0. This yields the infinite system of equations 



(n + a — e)b n - iX R ^fna n _i + X D \ / n+ la n+1 = 

(n + /3- e)a n + iX R \/n + lb n+1 + A D v^n-i = (5) 
for n > 0, with a_i = 0, 6_i = 0, and a = (1 + lu l /uj c )/2, (3 = (1 — uo L /uo c )/2. 

A. Case in which either = 0, or Ad = 0. 

When only the Rashba or Dresselhaus terms are considered, Eqs. (5) can be exactly 
solved. 14 > 30 ~ 32 For the sake of completeness, we give here the corresponding results. In the 
\d = case, combining Eqs. (5) one obtains 



[n + a — e)(n — 1 + f3 — e) — n X 







(n + a — e)(n — 1 + (3 — e) — n \ 2 R a n _i = , (6) 
either of which yields the energies 



e ± n =n± x \(l + ^-)\-Xln . (7) 



,4V UJr.J UJr 

One also obtains 

(n-l + P-4)a^_ 1 = -iX Ry /Eb^ , (8) 



which together with the normalization condition |a n 1i| + \bff\ = 1 solves exactly the 
problem [for n — 0, a_i = 0, bo = 1, and eo = |(1 + wl/w c )]. 

Eqs. (6) indicate that in the series expansion of the spinor only one and one foj 
coefficient appears. Specifically, 

In the limit of zero spin-orbit, the spinors \n d ) and \n u ) become and |^)^o^, 

respectively. The exact expressions for the Oj and bi coefficients entering Eq. (9) are easy to 
work out. Expressions valid up to X 2 RD order are given in the next subsection. 
The X R = case can be worked out similarly. One obtains the secular equation 

{n + (3-e)(n-l + a-e)-n~X 2 D = (10) 
4 



which yields 



One also obtains 



£ ± = n± x l- A [l- — ) +-Xln . 



{n-l + a-e^)b e 1 c_ 1 = -X D VEa e n n , 



J n-1\ 



(11) 

(12) 

1 solves exactly the 



which together with the normalization condition |a^ n | 2 
problem (in this case, for n — 0, = and ao = 1). 

Again, in the series expansion of the spinor \<f>), only one and one 6j coefficient appears: 



" +1 in+l) 



*n+l 



n„> = 



V 6;" +1 |n) 
and the same comments as before apply. 




(13) 



B. General case when \r ^ and Ad ^ 0. 



If both terms are simultaneously considered, the SO interaction couples the states of all 
Landau levels, and an exact analytical solution to Eqs. (5) is unknown, and likely does not 
exist. We are going to find an approximate solution that in the X RD /u c <C 1 limit coincides 
with the results of second order perturbation theory, i.e., it is valid up to \ 2 RD order, and 
it is quite accurate as compared with exact results obtained numerically. Combining Eqs. 
(5), one can write 

1 



T2 n 

n + a — e — A R 

R n- l+P-e 



X 2 



n 



n+ 1+fi-e 



b n = 



and 



-iXrXd 



n 



n(n-l) v /(n + l)(n + 2) 

0n-2 ; ; — — — — o n+ 2 



1 + 13 



n + 1 +0 - e 



n + f3 



A 



n + 1 



R 



n + 1 + a — £ 



I 2 



n 



n 



1 + a 



(14) 



-iX R X D 



n(n — 1) 



n — 1 + a — e 



a n -2 — 



(n+ l)(n + 2) 
n + 1 + a — e 



' O-n+2 



(15) 



The approximate solution is obtained by taking a„_ 2 = a n+2 = = b n+2 = in the above 
equations. This means that for each level |n), the SO interaction is allowed to couple it only 
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with the \n — 1) and \n + 1) levels. This solution, which consists of a \na) and a \n u ) spinor, 
is therefore obtained by solving first the secular, cubic equation 

(n+a-e)(n-l+P-e)(n+l+P-e) = \%n(n+l+P-e) + \' 2 D (n+l)(n-l+P-e) . (16) 

Together with the equations 

(n - 1 + (5 - £)a n _i = -i\ R y/nb n 
(n + l + fi- e)a n+ i = -\ D y/n + 1 b n (17) 

and the normalization condition |a n „i| 2 + |a n+ i| 2 + |6 n | 2 = 1, they determine the \rid) solution. 
The solution corresponding to the \n u ) spinor is obtained by solving the secular equation 

(n+(3-e)(n-l + a-s)(n+l+a-e) = X R (n+l)(n-l + a-e) + X 2 D n(n + l + a-e) . (18) 

Together with the equations 

(n — 1 + a — e)b n -i = —XoVncin 
(n + l + a- e)b n+1 = i\ R ^/n+ 1 a n (19) 

and \a n \ 2 + |6 n -i| 2 + |&„+i| 2 = 1, they determine the \na) solution. 

Since all the estimates available in the literature (see for example Refs. 15,21,22 and 
Refs. therein) yield \ 2 R D values of the order of 10 /xeV, and oo c in GaAs is of the order 
of the meV even at small Z?(~ 1 T), it is worth to examine the above solutions in the 
\ 2 RD = 2X 2 rd /uj c ^ 1 limit, in which the secular equations have solutions easy to interpret. 

To order \ 2 RD , the relevant solution to Eq. (16) containing both SO terms is 

e d = n + a + 2n—^ 2(n + 1) X ° , (20) 

LU C + UJ L LO c - UJ L 

that corresponds to the spinor \rid) 



, a e JLi \n — 1) + a^+i \n + 1) . 
\n d )= n+U (21) 



b £ n n \n) 



with coefficients 

UJ r . 



a^ n _! = i\ R Vn 



U c + 10 L 

a n n + i = -XdVtiTT — 

u c - UJ L 
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In the following, we will refer to this solution as to the quasi spin-down (qdown) solution, 
since in the zero spin-orbit coupling limit \rid) becomes |n)^ ]_ ^. Analogously, Eq. (18) has 
the solution 

e u n = w + /3-2(w + l) Xr + 2n X ° , (23) 

U) c + U) L U) C -U) L 



that corresponds to the spinor \n u ) 



I E U I \ \ 

^"-i|n-l>+£"+i|n + l> / 



(24) 



with coefficients 



K-i = *DVn 

uj c - uj l 

6„+i = iXnVn + 1 



n+1 ■ UJc + UJ L 

* = 1 -^» +1 )fe) 1 -^K^) a • (25) 

In the following, we will refer to this solution as to the quasi spin-up (qup) solution, since in 
the zero spin-orbit coupling limit \n u ) becomes \n) ^ q ^ . When either \ R or \ D are zero, Eqs. 
(21) and (24) reduce to the exact Eqs. (9) and (13), respectively, and the corresponding 
and bi coefficients, valid up to order A^d, can be extracted from Eqs. (22) and (25). These 
Eqs. show that a n and b n are of order 0(1), whereas a n ±\ and b n ±\ are of order 0(Ar 5 £>), 
and a n ±2 and b n ± 2 are of order 0{\ 2 R D ). This shows that the neglected terms in Eqs. (14) 
and (15) are of order 0(X RD ). 

The sp energies obtained from Eqs. (20) and (23), valid in the \\ D /uJc <^ 1 limit, are 

1\ , , a 2 UJ C n/ , lN , 2 



= (n + -H + + 2nA^ — - 2(n + l)\ 2 D 



2/ <- ■ ■ -Jt — \- - ■ / " - U 

2 uj c + ujl u c — &l 

El = (n + \)u c -f-2{n+ 1)A^^— + 2nA^^— . (26) 

Together with the structure of the associated spinors, Eqs. (21) and (24), this sp energy 
spectrum constitutes one of the main results of our work. By suitable differences of these 
energies, one may obtain the sp transition energies discussed in the next Sec. 

The above sp energies coincide with the ones that can be derived from second order 
perturbation theory with the standard expression 

p (2) _ I \- IH ~ i\RUJc(a + a- - a~o + ) + \ D uj c (a + a + + a~a_)|n)| 2 

^« "4 ^ E° - E° ' [ } 



where \n) = \n, |), \n, I) are the spin-up and spin-down eigenstates of the sp Hamiltonian 
\{a + a~ + a~a + )u) c - \uj l o z with eigenvalues -E°(T) — (n+ \)uj c - \ou L , and E°(l) = (n + 
\)u) c + respectively. 

The approximate solution Eq. (26) is very accurate in the high B limit (see below). It 
also carries an interesting information in the opposite limit of vanishing B. In this limit 
(cul,cu c <C Ajj D ), Eqs. (16) and (18) yield the solutions 

E d n = pu c {n\l+{n+l)\l] 

El = ^2uo c [{n + 1)A| + n\l\ (28) 

which show that, at B ~ 0, to order \ 2 R D the Landau levels are not split due to the SO inter- 
action, as one might naively infer from Eqs. (26). Another merit of the approximate solution 
is that it displays in a transparent way the interplay between the three spin-dependent inter- 
actions, namely Zeeman, Rashba and Dresselhaus. Such interplay has been also discussed in 
Ref. 22, in relation with the violation of the Larmor theorem due to the SO couplings, and 
in Ref. 33, where the Zeeman and SO interplay is discussed using the unitarily transformed 
Hamiltonian technique. Note also that in GaAs quantum wells, which are the object of appli- 
cation in this paper, due to the sign of g*, the lowest energy level is the qup one at the energy 
Eq = \uj c — \ujl — 2\ 2 r uj c / (uc + ul), containing the Rashba contribution alone, whereas the 
following level is the qdown one at the energy Eq = |u; c + — 2A|, u c / {uj c — u L ), containing 
the Dresselhaus contribution alone. For all the other levels both SO terms contribute to the 
level energies. 

We have assessed the quality of the above analytical solutions, Eqs. (26), by comparing 
them with exact numerical results for some particular cases. Indeed, the exact solution to 
Eqs. (5) can be obtained in the truncated space spanned by the lower M oscillator levels. 
Mathematically, Eqs. (5) are then cast into a linear eigenvalue problem of the type 

M(b)=e(b) , (29) 

where M is a 2jV x 2Af matrix, while a and b are column vectors made with the sets 
of coefficients {a n ,n = 0, . . . ,N — 1} and {b n ,n = 0, . . . ,Af — 1}, respectively. We have 
diagonalized M using a large enough Af to ensure good convergence in the lower eigenvalues. 
Fig. 1 displays a comparison of numerical (symbols) and analytical (solid lines) energies as 
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a function of the Rashba SO strength, for a fixed Dresselhaus strength, both in units of u c , 
namely y R = \ R /oo c and y B = \ 2 D /u) c = 0.01. The chosen values for yo and y R are within 
the expected range for a GaAs quantum well. For instance, if m\ RD /h 2 ~ 10/xeV and 
B ~ IT, (m\ RD /n 2 )/{hu c ) ~ 10- 2 . There is an excellent agreement between analytical and 
numerical results, differences starting to be visible only for strong Rashba intensities and 
high Landau bands. Actually, in Fig. 1 the largest value of the adimensional ratio between 
Rashba SO and cyclotron energy y R = \ R /oo c is 0.05, small enough to validate the analytical 
expression. Notice, however, that for larger y R values -not shown in the figure- i.e., for small 
enough B, Eqs. (26) no longer reproduces the numerical results. For GaAs this happens for 
magnetic fields below 0.1 T. Similarly, Fig. 2 displays a comparison of numerical (symbols) 
and analytical (solid lines) energies as a function of the Dresselhaus SO strength y^, for 
a fixed Rashba strength y R = \ 2 d /uj c = 0.01. For every Landau level, both figures show a 
crossing between the \n u ) state, which is at lower energy for y RjD <C 0.01 because g* < 0, 
and the \rid) state, that eventually lies lower in energy. This crossing is due to the interplay 
between both SO terms. 

III. SINGLE-PARTICLE LEVEL TRANSITIONS INDUCED BY APPLIED ELEC- 
TROMAGNETIC FIELDS. 

We can use the preceding results to study the sp transitions induced in the system 
by the interaction with a left-circular polarized electromagnetic wave propagating along 
the ^-direction, i.e., perpendicular to the plane of motion of the electrons, whose vector 
potential is A(t) = 2A(cos6i + sin #7), with 6 = out — qz. The sp interaction Hamiltonian 
J • A/c + g*^B s • (V x A), where J = ev/^/t, reads 

hi„t = (y-e i9 + v + e- w ) + ^g*fi B qA (a_e w + a + e- w ) , (30) 

where the velocity operator v± is defined as v± = — i[x ± iy, H] = P ± ± i\ R a± + \do' t . 
The Hamiltonian h int can be rewritten as 

h^t = ^r e A ^c (a-e w + a + e~ ie ) + ^g*fi B qA (W 9 + a + e~ ie ) , (31) 

where the operators a + and a~ acting on the spinor \<p) are 



a+ 


iX R 




a 


Ad 






, a = 






Id 


a+ 




-iX R 


a - 



(32) 
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In the dipole approximation (q rs 0), the charge-density excitation operator is v±. We note 
that, even in the presence of e-e interactions, this operator satisfies the f-sum rule: 

52(0\x=fiy\n)(n\-iv±\0) ='£uj nO \(n\x±iy\0)\ 2 = ho\[x=fiy,[H,x±iy]]\0) = 2N , (33) 

n n ^ 

where N is the electron number and uj ti o are the excitation energies. 

We consider next several useful examples of sp matrix elements involving the operators 
a + , which is proportional to v+, and cr_, and the qup and qdown sp states of Eqs. (21) and 
(24). For the operator a + , we can write in general 

= + i\ R Vi<h + W 2 *0i + ^*« + 02 , (34) 

and have to distinguish between qup-qup, qdown-qdown, qup-qdown, and qdown-qup tran- 
sitions. The qup-qup and qdown-qdown transitions represent the usual CR, and the qup- 
qdown and qdown-qup are related to spin-flip transitions. 

Let us start with the qup-qup and qdown-qdown transitions. To the order \ 2 RD they 
are dominated by the transition n — > n + 1 at the energies E% +1 — E% and E% +1 — E% with 



matrix elements |((n+ l) M |a + |n M )| = |((n + l)d\a + \nd) \ = \Jn + 1. The energy splitting of 
the cyclotron resonance is 



AE, 



CR 



2 ,, 2 U c 



4A^— 4A D 



(35) 

L) c + UJ L U) c - U) L 

The a + excitation operator also induces a qup-qdown transition with energy E% — E™ and 
matrix element | (nd\a + \n u ) \ = \dujl/(^c — ^>l)- This is a spin-flip transition. In particular, 
when n = it is related to the Larmor resonance at the energy 22 

AE L = ujl + 2 (\ 2 R — ^ \ 2 D — — — \ . (36) 

Note that the transition matrix element is linear in X D , and that in the presence of the 
Rashba interaction alone, a + causes no spin-flip transition. 

For the operator er_ one gets (ip\a-\<f)) = 2-^201 • The dominant transition is the spin- 
flip excitation at energy E% — E% with matrix element |(rid|<7_|n u )| = 2. The qup-qup and 
qdown-qdown cyclotron resonances at energies E% +1 — E* and E% +1 — E% are also excited 
with strengths \{(n + l) u |<7_|n u )| = \{(n+ l) d \o-\n d )\ = 2\ D y/n + 1 uj c / (uj c - uj l ). 

Other excitations that deserve some attention are those induced by the operators a + a± 
and a + a z . They are detected in inelastic light scattering experiments as spin dipole 
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resonances. 27 The operator a + a z excites the same cyclotron states as a + , at the ener- 
gies E% +1 — E% and E% +1 — E%, and with the same transition matrix element \/n + 1. 
In contrast, the operator a + a + mainly induces the transition from qdown to qup states 
at the energy E™ +1 — E%, whereas the operator a + <j- induces the transition from qup 
to qdown states at the energy E% +1 — E%- The transition matrix elements are given by 
|((n+ l) u \a + a + \nd) \ = \{(n + ^)dW + °'-\ n u) I = %V n + 1- We thus see that the dipole tran- 
sitions between Landau levels \n) and \n + 1) at 'unperturbed' energies E n+1 — E n are split 
by the SO interaction, an effect that under some circumstances may be observed, as will be 
discussed in Sec. V. 

IV. ELECTRON-ELECTRON INTERACTION AND SUM RULES 

In this section our aim is to discuss the role played by the e-e interaction in the physical 
processes in which SO effects can be important and, as a consequence, have a chance to be 
experimentally detected. Since we have obtained a spinor basis that includes the SO effects 
(SO basis), one might use it to diagonalize the Coulomb interaction. This has been done, 
for example, in Ref. 16, where the spinor basis Eq. (9) has been used to study the influence 
of the Rashba interaction on the incompressible Laughlin state. One could also use the 
SO basis Eqs. (21) and (24) to solve the random-phase-approximation equations, 34 or to 
study SO effects on the collective states of the quantum well in the adiabatic time-dependent 
local-spin-current density approximation. 22 ' 35 We have chosen a different way to incorporate 
interaction effects that, while being more approximate, it is accurate enough and allows one 
to obtain simple analytical expressions for the quantities of interest here. It is the sum rule 
approach, which is well suited to address the interplay between SO coupling and the e-e 
interaction in some relevant excitation processes. 

Let us firstly recall that in the absence of the SO coupling, two important theorems hold 
for the quantum well Hamiltonian H, in which the e-e interaction is included. They are the 
Kohn theorem 

3 3 

which tells us that, in photoabsorption experiments on quantum wells, a narrow absorption 
peak must appear at the cyclotron frequency ui = u c excited by the cyclotron operator 



11 



J2j Pj, an d the Larmor theorem 

lH,S-]=u L S- , (38) 

which states that in inelastic light scattering experiments at small transferred momentum, 
or in electron-spin resonance experiments, a narrow collective state must be excited by the 
Larmor operator S- = J2j at the Larmor frequency. These two modes are not influenced 
by the e-e interaction. Things radically change if we include in H the SO interaction. We 
then obtain 

[H, £ P+] =uj c J2 {P + + iA*<7+ + \ D a_)^ (39) 

3 3 

and 

[H, S-] = uj l S- + £ (2i\ R p-a g + 2X D P + a z )^ , (40) 

3 

which show that the SO interaction couples the cyclotron (dipole and spin dipole) and Lar- 
mor modes. We have studied in Sec III the effects of the SO coupling on these excitations 
in the absence of the Coulomb interaction. Now, we want to determine whether the pres- 
ence of both the SO and Coulomb interactions has an effect on the Larmor and cyclotron 
frequencies or, on the contrary, the results of the previous Sec. still hold. With this goal in 
mind, we introduce the following mixed sum rules 36 ' 37 

™t = ^E^o((O|^|0n)(0n|G t |O)±(O|Gt|0 n )(0 n |F|O)) 

" n 

= l - ((0\F(H - E o ) k tf\0) ± (0\G\H - E ) k F\0)) , (41) 

where |0) and \<p n ) are the exact gs and excited states of the full Hamiltonian H (including 
e-e interactions), and u n o = E n — E are the corresponding excitation energies. For k = — 3 
we obtain 

m " = \(0\[F,G^]\0) 
mt = i(0|[F,[if,Gt]]|0) 
m 2 - = l(0\[[F,H],[H,&}]\0) 

mj = ±(0\[[F,m,[H,[H,&\]]\0) . (42) 

Clearly, the more sum rules are known, the better knowledge of the Hamiltonian spectrum. 
With the four sum rules of Eq. (42) we can obtain information only on two excited states 
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-see below. Consequently, we will limit the analysis to the cases in which either the Rashba 
or Dresselhaus SO terms are present because, as one can see from Eq. (39) and Eq. (40) as 
well, in this case only two states would then be coupled by the corresponding SO interaction. 

Let us first consider the case where F = G = J2iP t r , i- e -j is the cyclotron operator. 
Evaluating the commutators in Eqs. (42) we have, to order X RD , 



m r . 



lUr, 



ran 



2Nu c 
2Nul 



2Nul 



2Nut 



1 - —(A 



(43) 



To obtain these Eqs. we have used that J2iP-T\fy = an d have assumed that the gs of the 
system is fully polarized, i.e., (0| X)i°^|0) — N. As such, these expressions are useless unless 
the left-hand side can be directly evaluated from the definition Eq. (41), and this evaluation 
yields a closed expression for the excitation energies and transition matrix elements. This 
is the case if we consider either of the Xr^d terms alone, because only two states are excited 
by the cyclotron operator = J2iPt acting on the gs |0). Dropping e.g. the Xr term, a 
straightforward calculation yields 

2uj r 



iii = 2Nu c 
2N— 



(UJr. 



\ 2 



2ul 



7T 2 



(LU, 



X 2 



LU 10 = LO c + 



2uj r 



- A D 

2co c x2 

D 



^20 = u L : — \" D , (44) 

where tti and tt 2 are the transition strengths to the cyclotron \<f> ni ) and Larmor \(f) n2 ) states, 
tti = I (0m I Si-Pj + |0)| 2 and 7r 2 = \{<f> n2 \ J2i -Pi + |0) | 2 , and u w , u 2 o are the respective excitation 
energies. This is in full agreement with the results of Sec. Ill, and shows that the e- 
e interaction does not affect the frequency and transition strengths of the cyclotron and 
Larmor resonances. 

The case A# = can be worked out similarly, and the same conclusion may be extracted. 
We recall and stress again the results obtained in the previous section, namely that when 
\d = 0, the Larmor state \<f> n2 ) is not excited by the cyclotron operator J2iPt ( n 2 turns 
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out to be zero). Alternatively, all previous calculations could have been carried out using for 
the Larmor operator, namely, F = G = J2i a +- Assuming again that (0| X^cr^O) = N, 

we obtain the same results and draw the same conclusions as before. This is a consequence 

of the structure of Eqs. (39) and (40). 

Using more sum rules, e.g. and mf, one may obtain information on other states that 

can be excited by the cyclotron operator J2iPt- Their consideration shows that the e-e 

interaction does not affect, to order neither the cyclotron nor the Larmor state, whose 

frequencies are the same as determined in Sec. Ill when both the SO terms are included in 

H. 

When the gs of the system has both qup and qdown occupied states, 21 the spin dipole 
operator J2i Pt a% z entering Eq. (40) excites a state at an energy u; c (l+/C) -see below, instead 
of oj c as it corresponds to the cyclotron (charge dipole) operator J2iPt> an d the results in 
Eq. (44) must be corrected for. This effect is not related to the SO interaction, and appears 
even in the absence of it. The spin dipole operator does not commute with the e-e interaction 
as the cyclotron operator J2i Pt does, and K, is precisely the contribution to the spin dipole 
operator m~l sum rule arising from the e-e interaction when one takes F — G — J2i Pf^z'- 

m+ = E^o|(0nlE^|O)| 2 

n i 

= ^(0|E^[^E^]]|0> = ^ C 2 (1 + /C) , (45) 

where 

^=2]^2<0|EV^(r y )(^-a2) 2 |0> . (46) 

K, can be extracted from inelastic light scattering experiments. 27 It turns out to be zero 
for fully polarized ground states, and small and negative -of the order of 10~ 2 - otherwise. 
Similarly, the spin flip dipole operators J2i Pt a ±i whose excitations can be also measured by 
inelastic light scattering, do not commute with the e-e interaction, which give rise to some 
energy corrections. It turns out that these corrections are equal for the three spin dipole 
operators J2i Pt a l,± because the value of K is the same for all them. Hence, the energy 
splittings among these excitations are not influenced by the e-e interaction, depending only 
on the Zeeman and SO energies as found and discussed at the end of Sec III. 

Finally, we want to comment on the consequences of the failure of the Kohn theorem due 
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to the SO coupling using the mf sum rule for F = J2iPi &\ an d G = J2iPi '■ 
mt = 5>no(0| £ Pr a l\ n }(n\J2P+\0} 

n i j 

= ^(o|[E^,[^E^ + ]]|o) = ^(oiE^lo) • (47) 

i j i 

This sum rule allows to study the interplay between charge and spin modes. If we cast it 
into a sum over 'spin dipole states' \m) and another over 'charge dipole states' \p), we obtain 



mi 



= E^o(oie^:ip)(pIE^ + io) 

P i j 

+ E^o(o|E^XHHE^ + lo) • (48) 



3 



If there is no SO coupling, Kohn's theorem holds, implying that (m| J2i?i |0) = 0. Thus, 
when the spin gs 2S Z = (0| X^cr^O) is not zero [otherwise, mi = from Eq. (47)], only the 
density modes would contribute to mi through the p-sum in Eq. (48). On the contrary, if 
the SO coupling is taken into account, Kohn's theorem is violated and the spin and charge 
dipole states are coupled to order X 2 R D , with (m| J2i -^ + |0) being now different from zero. 

To be more quantitative, let us assume that only one charge dipole state, the cyclotron 
state \p) at energy E\ = uj c + 0(\ 2 RD ), contributes to the first sum in Eq. (48), and only 
one spin dipole state \m), at energy E 2 = oj c (\ + /C) + 0(\ RD ), contributes to the second 
sum, where we have indicated by 0(X 2 RD ) the SO correction to the cyclotron and spin dipole 
energies. Let us define the mixed strengths 

ti = (o|E^XIp)(pIE^ + Io) 

i j 

n 2 = (o|E^XMHE^ + l > • (49) 

i j 

Evaluating the sum rules m$ and mi for the operators G = J2i Pf an d F = J2i Pf a l: one 
easily obtains 

n 2 = Wc(0 |l>j|0)|^ = u, c <0|I>i|0>^|^ . (50) 

Eqs. (50) explicitly show that if (0|X^Cz|0) — o, or if the SO coupling is neglected, the 
mixed strength tt 2 is zero, and the spin dipole state cannot be excited in photoabsorption 
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experiments. The strength 7r 2 is nonzero only at odd filling factors v [y = 2ir£ 2 n e , where 
is the magnetic length and n e is the electron density), for which 2S Z /N = 1/u. 
Besides, when the system is fully polarized at v — 1, the operators Y,iP t r and Y,i Pi a% z 
coincide and excite the same mode, so there is no splitting. The SO corrections 0(A 2 ) can 
be calculated by taking into account the occupation of the ground state, either using the 
sum rule approach of this Section, or the method of unitarily transforming the Hamiltonian, 
as described in Refs. 21,33. This calculation yields the energy splitting of the CR we discuss 
in the next Section. 



V. COMPARISON WITH EXPERIMENTS AND DISCUSSION 

An actual confrontation of the theoretical results we have obtained with the experiments 
is not an easy task because of the smallness of the SO effects, and because of the way they 
are presented in the available literature, which makes it extremely difficult to carry out a 
quantitative analysis of such a subtle effect. Thus, we have to satisfy ourselves with a semi- 
quantitative analysis, or to point out that these results are compatible with fairly rough 
estimated values of the SO coupling constants. We present now three such examples and a 
possible way to increase SO effects so that they could be easier to determine. 

Using unpolarized far-infrared radiation, Manger et al. 20 have measured the cyclotron 
resonance in GaAs quantum wells at different electron densities. The main finding of the 
experiment is a well resolved splitting of the CR for z/=3, 5, and 7, and no significant 
splitting for v = 1 and for even filling factors. We have seen that the SO interaction couples 
charge- density and spin-density excitations yielding the SO splitting of the CR given in Eq. 
(35). However, this expression, by itself, is unable to explain the filling factor dependence of 
the observed splitting, for which one has to bear in mind that the SO coupling between the 
J2i Pf and J2i Pi° l z operators is strongly enhanced when the spin gs is not zero, as explicitly 
shown in Eq. (50). We have also noted that /C contributes to the splitting. Eq. (35) has to 
be generalized to include these features. We obtain 



AE, 



CR 



4 (4^— - A 2 , — — — ) KU 



(51) 



where the factor 2S Z /N takes into account the actual sp contents of the gs. This equation, 
together with Eq. (50), embodies the theoretical explanation of the experimental findings. 20 
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In particular, it gives an appreciable splitting only for odd filling factors, for which the spin 
ground state S z is not zero. The analysis of the experimental splittings using the Eq. (51) 



recently used to reproduce the spin splitting in quantum dots 15 and wells. 22 This is, in our 
opinion, one of the most clear evidences of a crucial SO effect on a physical observable, 
because its absence would imply that the physical effect does not show up. 

The spin splitting of the first three Landau levels of a GaAs quantum well has been mea- 
sured in a magnetoresistivity experiment by Dobers et al. 24 We have shown in the previous 
sections that this splitting is not influenced by the e-e interaction, and that there is no spin 
splitting as B goes to zero [Eq. (28)]. Both facts are in agreement with the analysis of the 
experimental data, and with previous theoretical considerations 38 about the 5-dependence 
of the gyromagnetic factor g* , whose determination was the physical motivation of the mag- 
netoresistivity experiment presented in Ref. 24. These authors have derived a B- and n- 
dependent g* factor g*(B,n) = g$ — c{n + \)B, where g^ and c are fitting constants that 
depend on the actual quantum well. The possibility of a SO shift was not considered, and 
their chosen law for g* implies that the spin splitting energy AE n does depend on the Lan- 
dau level index n entering in a B 2 term, as they have AE n = \g* HbB\- A B dependence in 
g* is crucial to explain the experimental data, and also to reproduce them theoretically. 22 

For the spin splitting of the Landau levels we obtain 



-recall that ujl = \g*/dBB\- i.e., a splitting that increases with n because of the SO coupling. 
This SO correction has been worked out for the n = level in Ref. 22 using the equation 
of motion method. It is known that the experimental results 24 for n — 1 and 2 can be 
reproduced if g* depends on n and B, as already shown in that reference. We have verified 
that the n-dependence of g* cannot be mimicked by the n-dependence introduced by the SO 
interaction, Eq. (52). Recently, the analysis of g* has been extended to a wider magnetic 
field range using time-resolved Faraday rotation spectroscopy 28,29 

As a third example, we address the inelastic light scattering excitation of the spin dipole 
modes at v = 2 as measured in Ref. 27. For this filling factor, in the absence of SO coupling 
the spin-density inter-Landau level spectrum is expected to be a triplet mode 34 ' 39 excited 
by the three operators J2i -Pj +0 "z,± with energy splittings given by the Zeeman energy ojl- In 



yields values for the quantity m\\ 2 R — \ 2 D \/Ti 2 of about 30/xeV, in agreement with the ones 




(52) 
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the presence of SO interactions, we still expect a triplet mode to appear. Indeed, for v = 2 
we have S z — and the cyclotron and spin dipole modes excited by the operator J2i Pt a \ 
are decoupled, as previously discussed. Thus, for this operator only one single mode should 
be detected at an average energy oj = u> c (l + K). The other operators J2iPi i ~< J ±, yield the 
two other spin dipole modes at the energies 

?± rt , ir\ _i_ , . _l a ( \i LJc \2 ^ 



E ± = uj c (1 + JC)±lu l ±4(\ 2 r ^ X 2 D c —\ . (53) 

The splitting is thus symmetric and depends on the SO strengths. In the experiment, triplet 
excitations were observed in all measured samples up to electron densities corresponding to 
r s = 3.3 (we recall that r s = 1/ ^im e ). B was accordingly changed to keep the filling 
factor at v — 2. Only one triplet mode spectrum at B = 2.2 T was shown. From this 
spectrum, we infer that there is space for a ~ 5 — 10% SO effect on the splitting, assuming 
that at this fairly small magnetic field, g* is that of bulk GaAs, g* = —0.44. We have 
estimated that m\\ 2 R — \ 2 D \/fi 2 ~ lOfieV, in line with the previous findings. Systematic 
measurements, especially at high B where the splitting is larger, are called for to allow for 
a more quantitative analysis. 

Another unequivocal signature of spin-orbit effects in quantum wells would be the detec- 
tion of the Larmor state in photoabsorption experiments. The strength of this transition is 

A 2 

given by 2- n -( ^ ) 2 and only depends on the Dresselhaus SO coupling -see the comment 
immediately after Eq. (35). In most experiments, B is perpendicularly applied to the plane 
of motion of the electrons, and for GaAs the strength is so small that it has never been 
resolved. 

We finally discuss the effect of tilting the applied magnetic field using the expressions 
derived in the Appendix. Eq. (A. 3) can be used to obtain the splitting of the cyclotron 
resonance which generalizes Eq. (35) for tilted magnetic fields: 



AE, 



CR 



(C R V + C d Z) - ■ , I — - - (C R Z + C D V)- ! 



(54) 



l + \g*\m*S/2 v J l-\g*\m*S/2 

where C_r,d = m -^ijD/^ 2 > an d the tilting angle 9 enters the quantities V, Z, and S defined 
in the Appendix. Tilting effects might arise because of the 1 — \g*\m*S/2 denominator in 
the above equation, but sizeable effects on AE C r should only be expected for materials such 
that \g*\m* /2 is large. This is not the case for GaAs, but it is, e.g., for InAs and InSb, which 
have \g*\m*/2 = 0.169 and 0.355, respectively. For the latter case the dependence of AE C r 
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with the in-of-well field B x , with a fixed B z , is shown in Fig. 3. Notice that AE C r is sharply 
increased when B x exceeds a given value (IT for the parameters in Fig. 3), which is proving 
the strong enhancement of SO effects introduced by the horizontal component of the tilted 
field configuration. Figure 3 also shows the comparison with the exact diagonalization data 
(symbols), indicating that the analytical formula, Eq. (54), is accurate up to rather large 
tilting angles and for varying relative weights of Rashba and Dresselhaus terms. As a matter 
of fact, this analytical result does not depend on B z although, for the sake of comparison 
with the exact diagonalization, we have used B z = 1 T in Fig. 3. The evolution with B x is 
not always monotonous, especially for Cr > Cd where we find an initial decrease of AEcr 
with increasing B x , vanishing at B x ~ 0.8 T, and eventually increasing again. 
The tilting also affects the spin splitting of the Landau levels 

— = ^^+2(2n+l) L v + yD z) 1 _ ( z + yD V) i 

u c 2 1 + \g*\m*S/2 1 — \g* \m*S/2 

(55) 

which generalizes Eq. (52) for 6 ^ 0. As we have commented before, in a recent exper- 
iment where spin precession frequencies in a InGaAs quantum well have been measured 
using electrically detected electron-spin resonances, 29 a strong dependence of the effective 
gyromagnetic factor g e ^ on the applied tilted B has been found. In particular, at 6 — 45° 
g e H exhibits oscillations with B which indicate its sensitivity to the Landau level filling, and 
a coupling between spin and orbital eigenstates which is explicitly present in the spin-orbit 
term of Eq.(55). The effective g-factor that can be extracted from this equation at 6 — 45°, 
by taking the ratio 2AE n /(m*Suj c ), has the structure 



\g* f, (B,n)\ = \gZ\ + (n + ± 



(56) 



where the parametrization g* = g$ — C\{n + ^)B of Refs. 24,29 has been introduced in Eq. 
(55), and the C2 term is the SO contribution. For the smaller B values in the experiment, 
and for reasonable values of m\ 2 RD /h 2 ) of the order of 1-10/ieV, the SO contribution is 
important enough and should not be neglected; under these circumstances, time-resolved 
Faraday rotation spectroscopy could be sensible to Rashba and/or Dresselhaus spin-orbit 
effects. 
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VI. SUMMARY 



We have discussed the appearance of spin-orbit effects in magnetoresistivity and inelastic 
light scattering experiments on quantum wells. In particular, we have addressed SO effects 
on the splitting of the cyclotron resonance, on the sp Landau level spectrum, and on spin- 
density excitations. Our discussion has been based on the use of an analytical solution of 
the quantum well Hamiltonian valid up to second order in the SO coupling constants. The 
accuracy of this solution has been assessed comparing it with exact numerical diagonaliza- 
tions. 

We have carried out semi-quantitative comparisons with available experimental data, 
with the twofold aim of extracting the value of the SO coupling constants and of indicating 
possible manifestations of the SO interactions. We have also pointed out that tilting the 
-usually- perpendicularly applied magnetic field might enhance spin-orbit effects, making 
them easier to detect. 

APPENDIX 

In this Appendix we generalize some of the expressions derived in Sec. II to the case 
in which B has a in-of-well component, e.g., B = (B x , 0, B z ). The Zeeman term then 
becomes \g* //bB-<t= — \uJ z L {<J x t&a.6 + a z ), where we have introduced the zenithal angle 9, 
tan 6* = B x /B z , and the 'z-Larmor' frequency uj z l = \g*fisB z \, with uj z l /uj c = \g*\m*/2. The 
Schrodinger Eq. (4) then becomes 

\{a + a~ + aTa + ) — uj z l /(2uj c ) — e iX R a~ + X D a + — [uj z l /(2uj c )] tan# 

—iX R a + + X D a~ — [uo z L / (2uo c )} tan# |(a + a~ + a~a + ) + uo z L /(2uo c ) — e 

(A.l) 

The calculation proceeds as before, Eq. (5) becoming 

[n + a — e) b n — ^tan# a n — iXr^/ti a n _x + X^Vn + I a n+1 = 

a- (3 ~ 

(n + (3 - e) a n — tan 9 b n + iX R Vn + I b n+1 + X D ^/n 6 n _i = , (A. 2) 

where a — (1 + uo z L / uo c ) /2 and (3 — (1 - uo z L / uo c ) /2. 
The sp spectrum Eq. (26) becomes 

E d n ={n + ^ju c + ^S + 2n 



= 



U (X 2 R + X 2 D ) + (X 2 R V + X 2 D Z) 



UJr. 



LO, 



+ LUfS 
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-2(n+ 1) 



K 



-2(n + l) 



W (A^ + A^) + (A^ V + A 2 , Z) 



^ c - LU Z L S 



UJr. 



(A.3) 



w c + u z L S 

where we have defined S = l/cos0, W = sin 2 0/4, V = (l + cos#) 2 /4, and Z = (l-cos#) 2 /4. 
When 9 = 0, W = Z = 0, V = 1, and Eq. (A.3) reduces to Eq. (26). 
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FIG. 1: Lower energy levels for a GaAs quantum well as a function of the Rashba intensity 
VR = ^r/^c for a fixed Dresselhaus intensity yo = A^/a; c = 0.01. Solid lines are the analytical 
result, Eq. (26), while symbols correspond to the exact diagonalization, Eq. (29). 
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FIG. 2: Lower energy levels for a GaAs quantum well as a function of the Dresselhaus intensity 
Vd = ^b/wc f° r a fixed Rashba intensity ur = \ 2 r /oj c = 0.01. Solid lines are the analytical result, 
Eq. (26), while symbols correspond to the exact diagonalization, Eq. (29). 
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B x (T) 

FIG. 3: Splitting of the cyclotron resonance for an InSb quantum well (\g*\m*/2 = 0.355) as 
a function of the in-of-well field B x when B z = 1 T. Lines are the result from the analytical 
formula, Eq. (54), while symbols correspond to the exact diagonalization of Eq. (A.l). Defining 
Cr,d = m ^R d/^ 2 tne shown results are for: Cr = 30 /xeV and Cd = 10 ^eV, solid line and circles; 
Cr = 10 /xeV and Cd = 10 ^eV, long-dashed line and triangles; Cr = 10 /xeV and Cd = 30 /ueV, 
short-dashed line and squares. 
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